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1. Introduction 

Affine diffusions, introduced in the pioneering paper [9] by Duffie and Kan, are 
wideiy used in finance for modelling the term structure of interest rates. Their 
main attraction lies in the fact that they imply closed form expressions for bond 
prices. Affine diffusions are p-dimensional Markov processes that solve an affine 
stochastic differential equation (SDE) driven by a Brownian motion, i.e. an SDE 
with a drift n(x) and diffusion matrix 6{x), both affine functions in the argument 
x. There are three important issues in the theory of affine diffusions, to wit 

• stochastic invariance of a subset X of W, the state space, 

• the existence and uniqueness of strong solutions to the SDE with values 
in X, 

• the validity of the so-called affine transform formula for exponential mo- 
ments. 

Most of the theory that has recently been developed, concerns affine diffusions 
with a canonical state space R> x W~ m , henceforth referred to as the state 
space in standard canonical form, due to its tractable appearance which might 
ease the verification of possible technical conditions. The notion of a canonical 
state space has been introduced in [6] . Worth mentioning is the seminal paper [8] 
by Duffie, Filipovic and Schachermayer, who provide a complete characterization 
of (regular) affine processes, allowing jumps as well, under the assumption of a 
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standard canonical state space. Regarding the three issues mentioned above, for 
affine SDEs with canonical state space it is relatively easy to establish strong 
existence and uniqueness, as well as to derive conditions for invariance, see 
e.g. [6, 10]. Moreover, until recently, the afhne transform formula has only been 
fully verified for affine diffusions with a standard canonical state space, see [10]. 

The current paper together with a companion paper [19] contribute to the 
theory of affine diffusions with a non- canonical state space. We will characterize 
all affine diffusions, focussing our attention on polyhedrons (of which the stan- 
dard canonical state space is a special case) as well as quadratic state spaces 
(those of which the boundary is characterized by a quadratic function), though 
more is possible. For example, the matrix-valued affine processes and related 
Wishart processes treated in [4, 11] have the cone of positive semi-definite ma- 
trices as their state space. Our results extend the classification of [12] for the 
two-dimensional case to higher dimensions. In [12] it is shown that besides an in- 
tersection of halfspaces, also a parabolic state space is possible. We will see that 
in higher dimensions the quadratic state spaces are not limited to the parabolic 
ones; there exist also affine diffusions whose state space is a cone. 

The present paper concerns the first two of three mentioned issues for affine 
diffusions that live on a non-canonical state space. Results on the third one are 
presented in the companion paper [19], where we extend the results in [10] on 
the validity of the affine transform formula for canonical to general state spaces. 
Returning to the first issue, in the current paper we derive conditions for the drift 
and diffusion matrix on the boundary of X to ensure stochastic invariance for 
both the polyhedral and the quadratic state space. For the standard canonical 
state space these conditions are often called admissibility conditions, see [6] 
and [8]. The second issue, existence of a unique strong solution to an affine 
SDE, is in general not straightforward, as the square root of an affine matrix 
valued function 9 is not locally Lipschitz continuous for singular 9. This paper 
follows two approaches to solve this problem. 

The first one is by invoking a result by Yamada and Watanabc [21, Theo- 
rem 1], as is done in [6, 10]. This result is essentially only applicable for the 
standard canonical state space. Under invariance conditions though, we prove 
that a general polyhedral state space can be transformed in some kind of canon- 
ical form, not necessarily the standard one, for which the result by Yamada and 
Watanabe does apply. For a parabolic state space unique strong solutions can be 
similarly obtained by application of an appropriate modification of this result. 

The second way to obtain unique strong solutions is to impose conditions for 
invariance of [i £ P : 9{x) strictly positive definite} (also denoted by {9 > 0}), 
an approach followed in [9] for affine diffusions with a diagonalizable diffusion 
matrix and in [16] for matrix- valued diffusions. Strong existence and uniqueness 
is guaranteed, as the unique positive definite square root of 9 is locally Lipschitz 
continuous on {9 > 0}. In the present paper we derive invariance conditions for 
general state spaces, following the arguments in [16]. This enables us to obtain 
existence and uniqueness of affine diffusions whose state space is a cone. 

As a side note we mention that invariance of {9 > 0} is also important for 
applications. For example, in an affine term structure model one often desires 
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an affine structure of the underlying SDE under both the risk-neutral and the 
physical measure. For this purpose the invariance conditions for {6 > 0} are 
relevant in view of [19, Corollary A. 9], cf. [2, 7]. Therefore, we will provide these 
conditions not only for the cone but also for the polyhedral and parabolic state 
space. 

The remainder of this paper is organized as follows. After introducing in Sec- 
tion 2 some notation and defining affine SDEs and diffusions more carefully, as 
well as presenting a more detailed description of our aims, we discuss in Section 3 
stochastic invariance of the state space. Necessary conditions (admissibility con- 
ditions) on the behavior of the drift and the diffusion matrix on the boundary of 
a general closed convex state space are derived, whereas sufficient conditions are 
obtained for particular cases. In this section we also provide sufficient conditions 
for stochastic invariance of an open state space. 

The invariance conditions derived in Section 3 are used in Section 4 and Sec- 
tion 5 to characterize all affine diffusions with polyhedral respectively quadratic 
state space. For the former we also give sufficient conditions, extending those 
known from the literature ([3, 6, 10]), under which the diffusion matrix can be 
diagonalized. In particular we show that the classical model of [9] can be trans- 
formed into the canonical form of [6]. The results from convex analysis that we 
use in Section 4 are stated and proved in Appendix A. In Section 5 we show that 
for quadratic state spaces there are essentially only two types of state spaces 
possible, a (multidimensional) parabola and a cone. For each of these types we 
are able to give a full characterization of the possible diffusion matrices. 

2. Definitions, approach and notation 

Let p € N. We are given a p-dimcnsional stochastic differential equation 

dX t = n(X t )dt + a{X t )AW u (2.1) 

for continuous functions fx : R p — > M. p and a : W p — s- W xp that satisfy the linear 
growth condition 

\\fj.(x)\\ + \\<t(x)\\ < M(l + \\x\\), for all x E W, some M > 0. (2.2) 

By Theorems IV.2.3 and IV.2.4 in [13], for all initial conditions x G K p there 
exists a weak solution (A, W) to (2.1), that is, there exists a filtered probability 
space (0, T, (Jt),P) satisfying the usual conditions, with a p-dimcnsional F t - 
Brownian motion W and an adapted p-dimensional stochastic process A, such 
that Xq = xq a.s. and (2.1) holds. Let us also recall the result from Yamada 
and Watanabe [14, Theorem 21.14] that (2.1) has a unique strong solution if 
and only if weak existence and pathwise uniqueness holds. 
We use the following definitions. 

Definition 2.1. We call a measurable set X C W 1 stochastically invariant, if 
for all xo <G X, there exists a weak solution (A, W) to (2.1) with initial condition 
xq such that A t <E X almost surely, for all i > 0. 



P. Spreij and E. Veerman/Affine diffusions with non-canonical state space 



1 



Definition 2.2. The SDE (2.1) is called an affine SDE with state space X cW 
if 

1. it has a unique strong solution; 

2. X is stochastically invariant; 

3. the drift fi and diffusion matrix 9 = aa T are affine in x on X , i.e. 



for some a G W xp , b G W, symmetric JY 1 G W xp . 

The unique strong solution to an affine SDE is called an affine diffusion, which 
corresponds to the definition in [10, 19] in view of [19, Theorem 2.5]. 

Throughout the next sections we will address the following topics. Given a 
state space X : we determine all affine functions fi : W — > W and 9 : W xp — s- M p 
such that there exists a continuous square root a of 9 on X for which (2.1) is 
an affine SDE. To that end, the following three aspects have to be taken into 
consideration. 

• First, it is necessary that X C {9 > 0}, since 9(x) = a(x)a(x) T for x G X. 

• Second, for stochastic invariance, one has to impose conditions on /i and 
9 on the boundary of X , in order to prevent the solution X from leaving 
the state space X. 

• Third, one has to construct a square root a such that (2.1) admits a 
unique strong solution that stays in X. Remarkably, for the polyhedral and 
parabolic state space we consider, the conditions for stochastic invariance 
enable the construction of such a square root a, see Theorems 4.3 and 5.8. 
For the conical state space we obtain unique strong solutions by imposing 
conditions for invariance of {9 > 0}, see Theorem 5.16. 

Although obvious, it is worth noting that if X is an affine diffusion with drift 
yt{x), diffusion matrix 9(x) and state space X ', then LX + £ is an affine diffusion 
with drift Lfi(L~ 1 (x — £)), diffusion matrix L9(L~ 1 {x — l))L T and state space 
LX + i, for non-singular L G M pxp , I G M. p . Therefore, for the tasks as outlined 
above, it suffices to characterize all affine diffusions where the state space is in a 
certain "canonical" form (not to be confused with the standard canonical form), 
thereby obtaining all remaining diffusions by affine transformations. 

Remark 2.3. The distribution of X does not change with different choices of the 
square root a (as long as X stays in X for these choices) , since it is determined 
by the generator, which depends on a only through era 1 . Note however that 
if strong existence and uniqueness holds for one particular choice of a . it does 
not automatically hold for other choices. For instance take 9(x) = 1 in M and 
consider the 1-dimensional SDE dX t = <r(X t )dWt with aa 1 = 9. Existence and 
uniqueness of a strong solution holds when we take cr(x) = 1, while one only 
has a weak solution for the choice a(x) = sgn(x+), see [15, Example 5.3.5]. 



p 




(2.3) 
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Matrix notation 

The following notation regarding matrices and vectors is used throughout. Let 
p,q G N, P = {l,...,p}, Q = {l,...,q}, A G K px< ?, I C P, J C Q. Write 7 = 
{ii, ■ J = {ji, ■ ■ ■ ,j#j}, with i 1 < i 2 < . . . < i#i and ji < h < ■ ■ ■ < 

jff, j- Then Ajj denotes the (#/ x #J)-matrix with elements (Ajj)ki = Ai h j r If 
jfl = 1, say I = {i}, we write Aij instead. If J = Q then we write Aj instead. 
In particular, Ai denotes the i-th row of A. The j-th column is denoted by A? 
and the transpose of A is denoted by A T . The above notation is also used for 
matrix-valued functions <j>, e.g. (f> T (x) stands for (</>(cc)) T . 

For oi, . . . , a p G Kwe write diag(ai, . . . , a p ) for the p-dimensional diagonal 
matrix D with diagonal elements Da = ai, i G P. We also write diag(a) instead, 
where a denotes the vector with elements <Zj, sometimes explicitly denoted by 
a = vec(ai, . . . ,a p ). For a vector v we write \v\ for the vector with elements 
\vi\ and analogously y/\v\ denotes the vector with elements \Z|fj|. The identity 
matrix is written as I. We write [A] for the linear span of the row vectors of 
A. If A and B are two matrices with the same column dimension, then AJlB 
stands for L4]_L[B], i.e. AiBj = for all i and j. The unique positive semi- 
dchnite square root of a positive semi-definite matrix A is denoted by A 1 / 2 . 
For a square matrix A we write \A\ for (AA 1 ^) 1 / 2 . We will often make use of 
the fact that the matrix-valued function A H> l^l 1 / 2 is continuous, which is a 
consequence of [1, Theorem X.l.l]. 

3. Stochastic invariance 

In this section we obtain necessary and in some cases sufficient boundary con- 
ditions for stochastic invariance, see Definition 2.1. We first consider stochastic 
invariance of a closed convex set X C MP, for which we make use of the fact 
that it can be written as an intersection of halfspaces, i.e. 



with I some index set and U{ : W — > R : x h-> ^iX + Si, for some 7i G M lxp , 
Si G R. We denote the i-th boundary segment X n {ui = 0} of dX with dXi. 
The following proposition is partly proved in [10, Lemma B.l]. We give a more 
intuitive proof, involving an appropriate change of measure. 

Proposition 3.1. Let X C MP be a closed convex set given by (3.1) and assume 
X is stochastically invariant. Then necessarily it holds that 




(3.1) 



Vi e/,Vie dX t : ^a{x) = 
Vi G I,Vx G dXi : 7 ^(x) > 0. 



(3.2) 
(3.3) 



Proof. We give a proof by contradiction. Suppose there is an i G I and xo G dX\ 
such that (3.2) or (3.3) does not hold. Then there exists c G K such that 
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T) defined as 77(2) = Ji(J-{x) + cjia(x)a(x) T is negative in x . Let (X,W) 
be a weak solution to (2.1) with initial condition xq on some filtered prob- 
ability space (Q, J 7 , (J" t ),P). There exists a stopping time ti > such that 
L = £(c7iCr(X)l[o jri ] • W) is a martingale. Take T > arbitrarily, then we can 
change P into an equivalent probability measure Q on Ft by dQ = LrdF. By 
Girsanov's Theorem, W® defined by dV^ Q = dW t - ca(X t ) T jJl[ , T1 ] (t)dt, is 
a Brownian motion under Q on [0,T]. Hence X solves an SDE under Q for 
t G [0,T], namely 

dX t = (/i(X t ) + ca(X 4 ) ( x(X t ) T 7 . i T l [0 ^ l] (t))dt + a(X t )dW?. 

Let T2 > be a stopping time such that (c(X) ■ W®)tAr 2 is a Q-martingale and 
r;(X t ) < for t G [0, T2], Q-a.s. (this is possible since r)(Xo) = ??(^o) < and 
1 1-> Xt and 77 are continuous, P-a.s., whence Q-a.s.). Write r = ri A T2. Then it 
holds that r > 0, Q-a.s., and therefore 

/■TAt 

EQUj(XrAr) = + 7i E QX TAT = Ui(x ) + Eq / f7(X t )d£ < 0. 

Jo 

This implies Q(Vt £ [0,T] : Ui(X t ) > 0) < 1 and by equivalence of P and Q also 
P(Vt G [0,T] : Uj(Xt) > 0) < 1, which contradicts the stochastic invariance of 
X. □ 



If X C W given by (3.1) is a convex polyhedron, then / is finite, say I = 
{1, . . . , q}, for some q G N. We write 7 G R qxp for the matrix with row vectors 
7i, <5 G R g for the vector with elements Si and u : M. p — > R 9 : a; i-> 7.x + S. 
The necessary conditions obtained in Proposition 3.1 are sufficient when X is a 
convex polyhedron congruent to the canonical space R> x R p ~ m . 

Proposition 3.2. Let X C M. p be a convex polyhedron given by (3.1) with 
I — {1, ...,<?} /or some g G N cmc? suppose 7 /ias /itH row-rank. Then X is 
stochastically invariant if and only if (3.2) and (3.3) hold. 

Proof. This follows by [20, Remark 3.6], since X is congruent to the standard 
canonical state space, as 7 has full row-rank. □ 



Proposition 3.2 has the following corollary, which will turn out to be the 
building block for proving stochastic invariance for affine diffusions with general 
polyhedral state space, as considered in Section 4. 

Corollary 3.3. Let a G R pxp , b G W, c G E| and suppose the drift of (2.1) 
is given by 

H(x) = ax + 6, (3.4) 

and the diffusion coefficient by o~{x) = diag(y / ci|xi|, . . . , y/ Cp\x p \), for all x G 
R p . Then R> is stochastically invariant if and only if 



Mi G P,Mj G P\{i} : a tj > and b. t > 0. 



(3.5) 
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To obtain conditions for stochastic invariance of a quadratic state space, we 
note that if X C R p is a closed convex set with dX C = 0} for some C 1 - 
function $ : W -> R, then we can take (3.1) with I = dX, 7, = V$(i) and 
Si = — Tii, for i £ I, where we write V$ for the gradient of $ (written as a row 
vector). The necessary conditions for invariance (3.2) and (3.3) in this case read 

Vz G dX : V<f>{x)a(x) = 
Vx G <9A : V$(x)/i(a;) > 0. 

The next proposition gives necessary and sufficient conditions for a particular 
case. We write V 2 ^ = V T V<I> for the Hessian of a C 2 -function $. 

Proposition 3.4. Let $ : MP — > M : x i-> xi — /(xpun) /or some C 2 -function 
f : W- 1 -> R. TTien A" = {$ > 0} is invariant for (2.1) if and only if 

V$(x)ct(.t) = (3.6) 
V$(x)ju(a;) > -|tr (V 2 $(x)cr(x)o-(a;) T ), (3.7) 

/or all x £ dX = {<f> = 0}. 

Proof. Let (X, W) be a weak solution to (2.1) on some filtered probability space 
(fi, J 7 , (J-t),P) with initial condition xo € A. It holds that X t G A if and only if 
($(X t ),X P \ {1})t ) G R>o x RP- 1 . Ito's formula gives 

d$(X t ) = V$(X t )dX t + ±tr {V 2 $(X t )d(X t )) 

= (V$(X t )/i(X t ) + ^tr(yH(X t )a(X t )a(X t ) T ))dt 
+ V<f>(x)a(X t )dW t . 

Applying Proposition 3.2 gives the result. □ 



Note that if X is convex in Proposition 3.4, equivalcntly / is a convex func- 
tion, then V 2 / is positive semi-definite on its domain, whence V 2 $ is negative 
semi-definite and it follows that 

tr(V 2 <5>(X t )a(x)o-(x) T ) = tr {a(x) T V 2 ^(x)a(x)) < 0, for all x G M p . 

Thus condition (3.7) is stronger than (3.2), whence (3.2) and (3.3) are in general 
not sufficient for stochastic invariance. For more results we refer to [5] and [17]. 

We conclude this section by giving sufficient conditions for stochastic invari- 
ance of an open set. The idea behind the proof of the following result is taken 
from [16]. 

Proposition 3.5. Let <f> : R p — > R be a C 2 -function and suppose X° is a 
connected component (i.e. maximal connected subset) o/{$ > 0}. Then X° is 
invariant for (2.1) if there exists an open neighborhood O of dX such that 

V$(x)/j,(x) > -±tr(V 2 $(x)<9(x)) + i$(x)~ 1 V$(.T)6»(.x)(V$(.T)) T , (3.8) 

for all x in O D X° , where we write 6 = ercr T . 
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Proof. Let (X, W) be a weak solution to (2.1) on some filtered probability space 
(0, J 7 , (J-'t), P) with initial condition xo <G By Ito's formula it holds for 
t < t := inf{s > : $(X a ) = 0} that 

dlog$(X t ) = Stxy-MSOXt) - i$(X 4 )" 2 d($(X)) t 
= $(X t )- x (/(X t )df + a(X t )dW t ), 

where 

/(x) = V$(x)/z(x) + |(tr (V 2 $(x)6>(x)) - $(.t)~ 1 V$(x)6»(x)(V$(x)) t ). 

Suppose A := {to < oo} has positive probability. Since X does not explode, 
it holds for each uj G A that there exists e > such that Xt(uj) G On X° for 
toH - e < t < T (u;), whence f{X t {u)) > 0. Therefore, /„* $(X s )- 1 /(^s)ds 
does not tend to — oo on A for t f tq. This yields 

lim / $(X s )-V(X s )dW s = -oo, 

*tro Jo 

which is impossible, by the arguments of [16, Proposition 3.4]. □ 

Remark 3.6. Note that (3.8) implies (3.6) and is stronger than (3.7), as 6 = ao- T 
is positive semi-definite. 

Proposition 3.5 yields tractable conditions for stochastic invariance of an open 
set in case of an affine SDE. 

Proposition 3.7. Consider the situation of Proposition 3.5 and suppose ft and 
9 := u o~ are affine functions given by (2.3). Then X° is invariant for (2.1) if 

V$(x)6»(x) = $(x)w T , for some constant v <G W, (3.9) 
p 

V$(x)(/i(x) - i ^(A 1 ) 1 ) > 0, for all xeX, (3.10) 

where {A % ) % denotes the i-th column of the matrix A 1 . 
Proof. In view of Proposition 3.5 it suffices to show that 

p 

tr(V 2 $(x)6»(x)) - $(x)" 1 V$(x)6»(x)(V$(x)) T = -V$(x) ^(A*)*, 

under the assumption that V$(x)#(x) = $(x)w T , for some constant vector 
v G MP. Differentiating the right- respectively the left-hand side of the latter 
and applying the product rule yields 

V T $(x)w T = V T (V$(x)6»(x)) = V 2 $(x)6»(x) + V$(x)V T 6»(x), 
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where V$(x)V T 0(x) is short-hand notation for the matrix with row vectors 
d 

V$(x)—9(x). Hence 

oxi 

tr (V 2 $(a;)6>(») = V$(x)v - tr (V$(a;)V T <9(x)). 
The result follows, since V$(x)v = $(x)~ 1 \7 $(x)9(x)(\7 $(x)) T and 

P ( \ P 

tr (V$(x)V T 9(x)) = V$(a;)^-^ = V$(x) ^(A 4 ) 4 . 

i=l Xl i=l 

□ 

4. Polyhedral state space 
4-1- General diffusion matrix 

Throughout this section, we assume X is a polyhedron given by 

<? 

#=f|K>0}, (4.1) 

i=l 

for some affine function u : MM R p : i h> 7s + (5, 7 £ E« xp , <5 e R 9 , some 
9 £ N. We write Q = {1, . . . , g} and assume Q is minimal in the sense that 
P| ie g,{?ij > 0} 7^ A" for all Q' C Q with Q' ^ Q. In addition we assume /j, is of 
the form (3.4) and we are given an afhne function 9 by 

v 

9 :W ^W xp : x^ A Q + J2^x l: (4.2) 

i=l 

for some symmetric A 1 £ MP xp and we assume 

^ X° C X C {x e W : 6{x) > 0} =: V. (4.3) 

We will often make use of the fact that for a symmetric matrix S and vector v 
it holds that Sv — is equivalent with v Sv = 0. In particular if 9 = <t<t t on 
X, then (3.2) is equivalent to 

V« e/,Vie cU; : 7 i 0(x) = 0. 

In Section 3 we have derived necessary boundary conditions on ft and 9 to 
have X stochastically invariant. We show in this subsection that for a polyhedral 
state space, these conditions arc also sufficient for the existence of a square root 
a of 9 on X such that (2.1) is an afhne SDE. In the next proposition we prove 
that 9 can be transformed in a block-diagonal form, which we use in the proof 
of Theorem 4.3 to construct the square root a. The proposition improves upon 
the results given in the appendix of [9] and generalizes [10, Lemma 7.1]. 
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Proposition 4.1. Let X C W be a convex polyhedron given by (4-1) and sat- 
isfying (4-3) and let 9 be given by (4-2). Assume 

Vi G Q, Vx G dX t : ~/i6(x) = 0. (4.4) 

Then there exists a non-singular L G W xp and a vector £ G W such that 

/or some index sets M = {1, . . . , m} ; iV = {m+1, . . . , m+n} and affine function 
. In addition, we have 

LX + £ = K|' x C x E p - m -™, (4.6) 

for some convex polyhedron C = Di=i {^i — 0} C K.> wii/i = Xi for 

i < n. 

Proof. We divide the proof into a couple of steps. 

Step 1. There exists B eU qxp such that for i G Q it holds that 

7 ^(x) = P^x), (4.7) 
P j7 7 > 0, if Bj ^ (4.8) 
B a ] = 0, for j G (4.9) 

This is shown as follows. Fix i £ Q. By (4.4) and Proposition A. 4 there exists 
Bj G M lxp such that (x) = BijUi(x) for j G P. By assumption there exists 
xq E X° . It holds that iti(xo) > 0, so we can write 

P, = Ui (x o )S0(x o ). (4.10) 

By positive semi-definiteness of 0(xq), it holds that BijJ > 0. We have B^J = 
if and only if ~f t 8(x ) = 0, i.e. B t = 0. This yields (4.8). Moreover, if j G Q\{i}, 
then by symmetry of 8 it holds that 

B l jju l (x) = 7i0(x)'y] = Jj0{xhj = B^Ju^x). 

This implies B^J = 0, since Q is minimal. 
Step 2. It holds that 

rank 7j\/ = rankP^/ = m (4-11) 

K lxp = [7m] © N 1 , (4.12) 

where 

M = {i G g : Bi ^ 0}. (4.13) 

Indeed, by (4.9) we have 

B M1 T M = diag(P l7l T , • ■ • , B ml 1), (4.14) 
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which has full rank m by (4.8). This implies (4.11) as well as [7m] ^[Bm] 1 - = {0}, 
which yields (4.12). 
Step 3. It holds that 

[5m] 1 = [is] © [v], (4-15) 
[1q\m] = [is], (4.16) 

where N C Q\M is such that imun has full rank, r&nkjMuN = rank 7, for 
some 77 G ]R(p _ ™ 1_ ") x ?' j with n = #N. The first equality follows immediately 
from (4.11) and (4.12), while the second holds since [iq\m] C [Bm] 1 " by (4.9) 
and [7mujv] = [l]- 

Step 4- Permute indices such that M = {1, . . . , m}, N = {m + 1, . . . , m + n}. 
We define L G R pxp by 

Lmun = Imun (4-17) 
Lp\(MiiN) = V- (4-18) 

Then £ is non-singular, by (4.12) and (4.15). Moreover, the previous steps yield 

- t /diag(ciui(a;), . . . ,c m u m (x),0 N ) 



L0(x)L 1 



^ $(x) 



with Ci — BciJ and $ defined by $(x) = n9(x)ri T . It holds that Cj > by 
(4.8). By rescaling tij we may assume c, = 1. Then we take £ G K p such that 

S'tep 5. It remains to show that we can write &(x) as an affine function 
of Umus{x). This is an immediate consequence of the assumption that X C 
{9 > 0}, as this yields that 8(x) only depends on u(x), which is a function of 
u MuN (x). □ 

Remark 4.2. For the canonical state space as treated in [8] and [10] one has 
(4.6) with C = M| and 



*(xmus) =B q +J2 B ' x ^ 



with B l positive semi-definite. 

Theorem 4.3. Let X C K p be given by (4-1) and satisfying (4-3). There exists 
an affine SDE with drift fx, diffusion matrix 9 and polyhedral state space X if 
and only if 

Vi G Q,Vx G dXi : ii6(x) = 0, (4.19) 
\/i € Q,\/x € dXi : nn(x) > Q. (4.20) 

Proof. The "only if" -part is Proposition 3.1. We prove the "if" -part. 
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Assume (4.19) and (4.20). By Proposition 4.1 we may assume 

e{x) = ^( 7 M t(i 0j fo[x ^ (421) 

and 

X = M^ x C x R p -" 1 -™, (4.22) 

for some index sets M = {1, . . . , m}, N = {m + 1, . . . , m + n}, an affine function 
\P and some convex polyhedron C C K> , sucn that umun(x) — %mun- As a 
square root of 9 on X we take 

.W = Wx)|^ = (^(v«,o») |t(iif w)|1/2 ). (4 . 23) 

It remains to show that there exists a unique strong solution X to (2.1) which 
stays in X for all xq g A\ 

Let xo g K p be arbitrary. By continuity of the coefficients ft and a and satis- 
faction of the linear growth condition (2.2), there exists a weak solution (X, W) 
to (2.1) on some filtered probability space (£7, J 7 , (Ft), P). To show stochastic 
invariance of X, note that by condition (4.20) and Proposition A. 2 there ex- 
ist a g R 9 * 9 with a l} > for i,j g Q with i ^ j, and 6 g R?, such that 
7/i(x) = a«(i) + 6. Hence u(X t ) satisfies the g-dimcnsional SDE 

du(Xt) = (au(X t ) + b)dt + diag( v /|A A «|, QW )dW u (4.24) 

with initial condition u(xq), where W is a ^-dimensional Brownian motion with 
Wm = Wm (possibly after extending ft). For this SDE strong existence and 
uniqueness holds by [21, Theorem 1]. Together with Corollary 3.3 this yields 
u(X t ) g R q >0 , P-a.s., for all t > if u{x ) g E^ . In other words, X t g X, P-a.s., 
for all t > if xq € X. Thus X is stochastically invariant. 

We now show pathwise uniqueness for (2.1). Therefore, let (A, W) be an- 
other solution on the same filtered probability space (H, J 7 , (.F t ),P) with initial 
condition xq. Pathwise uniqueness for (4.24) implies that Xmun = umun(X) = 
umun(X) = Xmun a.s. Write R = P\(MU N). Since Xr does not appear in 
the diffusion part of the SDE, we have 

d(A_R jt — Xrj) = a_Rfl(Afl jt — A_R :t )dt, X,r : o — Xr$ = 0. 

So Xr — Xr solves a linear ODE which has as its unique solution, whence 
Xr = Xr a.s. and the result follows. □ 

By an inspection of the proof of Theorem 4.3 we see that if X solves an 
affine SDE with polyhedral state space X given by (4.6), then u(X) solves the 
affine SDE (4.24) with admissible parameters in the sense of [6] (that is, a has 
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non- negative off-diagonal elements and b € K> )- Extending the dimension by 
considering (u(X), Xr) instead of X, we get another affine SDE with state space 
u(X) x W (with r = and diffusion matrix 9 (say). Now if 9 can be written 
as an afhnc transformation of u with positive semi-definite matrices, then 9 is 
positive semi-definite on the whole of R> x W . In that case, u(X) x W can be 
enlarged to the canonical state space R> x R r and the resulting SDE is of the 
canonical form as in [10]. We elaborate on this in the next subsection. 

4-2. Diagonalizable diffusion matrix 

In [3] it is shown that the diffusion matrix 9(x) of an affine SDE with a standard 
canonical state space cannot be diagonalizcd in general. That is, there does 
not exist a non-singular matrix L such that L9(L~ 1 x)L r is diagonal. In this 
subsection we show that for a large class of polyhedral state spaces, including 
the standard canonical state space, diagonalization of the diffusion matrix is 
still possible in a different way, by extending the dimension. We also provide 
sufficient conditions for this as well as we give an example of an affine diffusion 
whose diffusion matrix is not diagonalizable. 

Proposition 4.4. Let X be an affine diffusion with drift \i, diffusion matrix 9 
and polyhedral state space X. Then X is in distribution equal to an affine trans- 
formation of an affine diffusion with diagonal diffusion matrix and canonical 
state space of the form R™ x W~ m if and only if 

q 

9{x) = B° + B l u,{x), for some positive semi-definite B % . (4-25) 
i=i 

Proof. We first prove the "only if'-part. Suppose X = LX + 1 for some matrix 
L and vector £, where X is an affine diffusion with diagonal diffusion matrix 
9 and polyhedral state space X. Then it holds that L9(x)L T = 6(Lx + £) and 
X = LX + £. Since 9(x) is diagonal, we have 

9{Lx + £) = L9(x)L T = ^ d 4 (x)L l (i 8 ) T , (4.26) 

i 

for some affine real-valued functions di- Note that L l (L l ) T > and di(x) > 
for x £ X. To show that 9 is of the form (4.25), it suffices to write di as 
di{x) = V , A, , n,\ L.v + 1) + a for non-negative A,j and c,;. Note that (4.26) 
yields dk{x + y) = di(x) for y £ kerL (otherwise, replace di{x) by di(H(x)), 
where II denotes the projection onto (kerL) T ). Hence the affine map 

<f>i : Lx + £ n- di (x) 

is well-defined and since 4>i{y) > for y € X = LX + £, Proposition A.l yields 
the result. 
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Next we prove the "if'-part. Suppose (4.25) holds. By Proposition 4.1 we 
may assume 9 is of the form (4.21) with 

1 

$(x M un) = A + A 'u t (x), 

for some positive semi-definite A 1 . In this case (u(X), Xp\ (mun)) IS an affinc 
diffusion where u(X) satisfies (4.24). Therefore, we assume without loss of gen- 
erality that u(x) = xmun, so that Q = MUN, q = m + n and X = R?, xW~ q . 
Since A* > 0, its unique positive semi-definite square root (A 1 ) 1 / 2 exists. We 
write 

A i/2 = ( (A 0)i/2 (A i)i/2 (A«)Va) 
w(xq) = vec(lQ,a:il(3, . . .,x q l Q ). 

Now note that 

Hxmun) = A 1 /2 d iag(«;( a;Q ))(A 1 /2)T j 
so X is an affine diffusion with drift /i and diffusion matrix 



9(x) 



I \ fdiag{x M ,0 N ) 
k 1 ' 2 ) \ diag(w(x Q )) J ^0 A 1 / 2 



We can diagonalize this by extending the dimension. Define a non-singular 
square matrix T by 

' A V2 r 

I 0; 



T = 



and let X be an affinc diffusion with drift (fj,(x),0) and diffusion matrix 

I 0\ fdmg{x M ,0 N ) \fl X 

T) \ &mg{w{x Q ),Q Q )) \Q T 

and with the law of Xp equal to the law of Xo. Then 

'I 



T- 1 ) X 

solves an affinc SDE with diagonal diffusion matrix diag(xA/, O^r, w(xq), 0q) 
and it is easy to check that Xp satisfies an affine SDE with drift fi, diffusion 
matrix 9 and initial condition the law of Xq. Hence Xp is in distribution equal 
to X in view of Remark (2.3), which yields the result. □ 



The relevance of diagonalizablc diffusion matrices 9 is elucidated in affine 
term structure models. In such models the short rate is an affine transforma- 
tion of an affine diffusion X, the state factor. In view of Proposition 4.4, for 
an unobservable state factor X, we may assume without loss of generality that 
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9 is diagonal when it is of the form (4.25). In particular this applies to afhnc 
diffusions in canonical form. This observation reveals that affinc diffusions with 
canonical state space and a non-diagonal diffusion matrix have the same po- 
tential as those with a diagonal diffusion matrix, which answers the implicit 
question in the concluding section of [3]. 

The following example shows that there exist affine diffusions with non- 
canonical polyhedral state space that do not meet the criteria of Proposition 4.4. 



Example 4.5. Consider the polyhedron X = C x I 2 C R , where we take 

3 
2 ' 



C = f) i=1 {Ui > 0} C R> , with m(x) = xi, u 2 {x) = x 2 , u 3 (x) = x\ + x 2 - 
Let fi and 9 be given as 



+ l \ „,W^„«, Qf~\ - ( dia g(°^) 

$( x ) 



Hn(x) = + }^ > ^p\n arbitrary , 9{x) 

with N = {1,2} and 



$(x) 



X! + | 1 
1 x 2 + I 



Then \x and 9 fulfil (4.19) and (4.20), so by Theorem 4.3 there exists an affine 
diffusion with state space X, drift \i and diffusion matrix 9. However, one can 
show that 9 is not of the form (4.25). 

We now give sufficient conditions for (4.25). In Proposition 4.6 below we prove 
that (4.25) does not only hold under full-row rankness of 7 (equivalent to the 
canonical state space) but also under the weaker condition of full row-rankness 
of (6 7) and an additional assumption. 

Proposition 4.6. Let X C R p be given by (4-1) and satisfying (4-3). Suppose 
that either 

(i) 7 has full row-rank or that 

(ii) (8 7) has full row-rank and for all i G Q, x € R p it holds that Uj(x) = 
for all j € Q\{i} implies Ui(x) > 0. 

Then (4.25) holds. 

Proof. The case when 7 has full row-rank is easy, since the state space can be 
transformed in the canonical form R™ x W~ m . We consider the case when 7 
has not full row-rank. First we assume 8% > for all i £ Q. 

Let T € r(p+ 1 ) x (p+ 1 ) be a non-singular matrix such that Tq = (8 7). 
We extend the dimension of 8 and 7 by writing T = (8 7) . We also write 
Ui(x) = jiX + Si, x €E R p , for all i < p + 1, in other words 



Ui{x) = Ti 

Let N = r _1 . Then we have 



xi ^ \x) 

7 v ' »=i 



= r _i r = V N l Ui (x), for x e 
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Let A = (A A 1 ... A p ) and abuse notation by writing 

P+i 

Ay = J^A i - 1 y i , for y G « p+1 . 

i=l 

Then we can write 

0(x)=A\J e ) =J2(AN*) Ul (x). 

Define A' = AN' 1 for i < p + 1. It is sufficient to prove that A 1 > for i < q and 
A' = for i > q. 

Let i < q and write N* = (N{, N$) with N{ 6 1, iV^ € W. We consider two 
cases. 

Case 1. Suppose N{ + 0. Then N l 2 /N{ G {(9 > 0} =: V. Indeed, for j g Q we 
have 

%W^) = ^r,^ = ] ^i {i=i} . 

So Uj(Ni/Ni) = for j g Q\{i}, whence by assumption Ui(Ni/N{) > 0, i.e. 
iV} > and iVl/iVJ g Af c V. It follows that 

< 6{Ni/N\) = —AN 1 = — A\ 
- \ 2/ v N * N i » 

and thus A* > 0, as N{ > 0. 

Case 2. Suppose N\ = 0. Then we have 

A* = AN< = ^eA^ /e )=ti T em/e). 

So we have to prove N\je g V for e small enough. For j g Q it holds that 

This gives 

Uj(Ni/e) = Sj + l{i=j} /e > 0, for j g Q, e > 0, 

since Sj > for j g Q by assumption. Hence N^/e G <Y C I? for e > and thus 
A* > 0. 

We have just shown that A* > for i < q. We now show that A 1 = for 
i > q. Let i > q and take xo g X. Write R = P U {p + 1}. Since rank7 = p 
and rank7 fl \{jj = p—1, there exists £ G R p such that 7j£ = for all j g 
and ^ 0. Then for all fc G R we have Uj(xo + fc£) = Uj(xo) for j G in 
particular for all j G Q. Hence xo + kt; £ X for all fc G R. Therefore, 

p+i 

< 0(x + k£) = ^2 A. j Uj(x + fc£) = 6(x ) + fcA 4 7^, for all fc G R. 

3=1 
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Since 7^ ^ 0, it follows that A' = 0. 

Now consider the general case without the restriction Si > for all i £ Q. 
Let xo £ X and for i £ Q define 

Ui : MP — > K : x \-¥ jiX + Si, 

by Ui(x) = Ui(x + Xo), i£ R p , i.e. Si = Si + jiXo and = 7$. Moreover, write 
9(x) = 6{x + x Q ), x £ W. Then Si = u,(0) = Ui(x ) > 0, i £ Q. Note that if 
uj(x) = for all j £ Q\{i}, then Uj(x + x ) = for all j £ Q\{i} and hence 
Ui(x) = Ui{x + xo) > by assumption. Moreover, 

(<* 7) = ( 5 + 7^o 7) = (<* 7) , 

which has full row-rank. Therefore we are in the previous situation, for which 
we have proved the existence of positive semi-definite A 1 £ R pxp for i £ Q U {0} 
such that 

<? 1 
9(x + x ) =9(x) = A° + ^A i u i (a;) = A + ^ A ?; w 4 (a; + &„). 

i=l i=l 

This gives the result. □ 

Corollary 4.7. Lei A" C M p &e given by (4-1)- Suppose X is contained in 
a polyhedron C C {9 > 0} which meets the criteria of Proposition 4-6- Then 
(4.25) holds. 

Proof. Write C = 01=1 { w « — ^ 0}> f° r some afiine functions Vi, with r G N. By 
Proposition 4.6 it holds that 

r 

0(a:) = B" +J2 Bi Vi(x), (4.27) 

i=l 

for some positive semi-definite B l . By Proposition A.l there exist Xij > 0, Cj > 
such that 

r 

Plugging this in into (4.27) we get the result. □ 

Example 4.8. In the 2-dimcnsional case, the polyhedrons which satisfy the 
conditions of Proposition 4.6 are the "triangles" (including those with vertices 
and edges in "infinity", like {0 < xi < 1} n {x 2 > 0} and {xi > 0}). Thus 
by Corollary 4.7, if X is contained in a triangle that is a subset of {9 > 0}, 
then (4.25) holds and X can be transformed in canonical form (see the proof 
of Proposition 4.4). However, this condition is sufficient, but not necessary. For 
example let 

6 W = ( X 1 xJ = (l o) + (o o) Xl + (o l)* 2 ' 
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and take u\(x) = 2x\ — x%, U2(%) — ~\ x i + x 2, ^3(0;) = — § + x\ + x-i- Then 
{0 > 0} = {x G K| : x 2 > ±} and X = ff i=1 {ui > 0} C {9 > 0}, but X is not 
contained in a triangle which is a subset of {9 > 0}. Still we can write 9 as an 
affine transformation of the Mi's with positive semi-definite coefficients, namely 

„ 1 A 1\ 1/4 2\ 1/2 4\ 1/2 -2\ 
9= 2{l lj + 9 U lJ Ml + H4 8j M2+ 9U 2j M3 ' 

^.5. Classical model 

In this subsection we revisit the classical model as introduced by Duffie and 
Kan in [9]. We assume 9 is of the form 

v 

6 = S diag(w)E T = ^ E^E*) 1 "^, (4.28) 

i=l 

with E G R pxp , u : R p -4 P : i 4 fix + a for some fi G R pxp , a G M p . In 
addition we assume X C HiLiI^ — 0} an d •'f 7^ 0. Under conditions (4.19) 
and (4.20), the proof of Theorem 4.3 constructs a square root a of 9 on X such 
that (2.1) has a unique strong solution. We show that the natural choice 

<7 = Ediag(^[) = (Si/M ■•■ £ p /W) 

also gives strong existence and uniqueness for (2.1). This is not immediately 
clear in view of Remark 2.3. 

Proposition 4.9. Suppose 9 is of the form (4-28) and assume X° / I, Z C 
P| p =1 {i>i > 0} and conditions (4-19) and (4-20) are met. Then (2.1) is an affine 
SDE fora = Y, diag( A /R) • 

Proof. Since conditions (4.19) and (4.20) are sufficient for invariance by Theo- 
rem 4.3, it suffices to prove existence and uniqueness of a strong solution. 

Note that X solves an affine SDE with a(x) = X diag(y|i; (x)|) if and only 
if LX + £ solves an affine SDE with cr(x) = LS diag(-y/|?7(x)|), where v(x) = 
v(L^ 1 (x — £)), for some non-singular matrix L and vector I. By Proposition 4.1 
we can therefore assume without loss of generality that 9 is of the form (4.21) 
and X of the form (4.22). Since all square roots of a positive semi-definite matrix 
are related by an orthogonal transformation, we have 

s <^^> - ( dia€(V f '°" ) (.(.JU*) <« 9 > 

for x € X, with 0(x) an orthogonal matrix, possibly depending on x. We show 
that there exists a matrix U G R( m + n ) x P with orthonormal rows such that 



Sa/uw diag(\/\v(x)\) = di&g(y/\x M \, N )U, for all x G W. 



P. Spreij and E. Veerman/Affine diffusions with non-canonical state space 



19 



Let xq £ X° and define U = Omun(xq). We have to show that 



\Vj(x)\ = s/WWij, for all i £ M, all j. (4.30) 

For i £ M we have 



S^y ^-(x) = \fx~iOij (x) , for all j and all x G A". 
If Ejj 7^ 0, then 

Uj(jc) = S^ 2 XiO?-(x), for all x € X, 

which yields Vj(x) = ex.; for some c > 0, by Proposition A. 4. If c > 0, then 
Oij(x) is constant, whence equal to Uij and (4.30) follows. If c = or E^- = 0, 
then Oij(x) = for Xi ^ 0, whence Uij = and again (4.30) holds. 

Let (X,W) be a weak solution to (2.1) on some filtered probability space 
(fi, J 7 , (J-t),P). Then [/W is an (m + n)-dimcnsional Brownian motion and it 
follows that Xmun solves an SDE with diffusion part diag(y^XM|, Oat). The 
strong existence and uniqueness for (2.1) follows along the same lines as the 
proof of Theorem 4.3, as o~(x) only depends on xmun by Proposition A.l. □ 



In the case that {9 > 0} 7^ 0, we can strengthen the result of Proposition 4.4. 

Proposition 4.10. Consider the situation of Proposition 4-9. If {9 > 0} 7^ 0, 

then the solution to (2.1) can be obtained by a bijective affine transformation of 
an affine diffusion with diagonal diffusion matrix and canonical state space of 
the form W£ a x RP~ m . 

Proof. As in the proof of Proposition 4.9 we may assume 9 is of the form (4.21) 
and X of the form (4.22). Since {9 > 0} ^ we have N = 0. Take xq £ X° = 
R> x W~ m with xo,i = 1 for i £ M . By strict positive definiteness of ^(xo.m) 
we can apply the linear transformation 

hlM \ 

(f(x„ iJ( ))-^J x ' 

so that we may assume without loss of generality that ^(xq^m) = h i-e. 9{xq) = 
I. Since v;(xo) > for all i, we can substitute Vi/vi(xo) for Vi and y / w^(xo)£ I 
for E 4 , which does not affect <r(x) and yields Vi(xo) = 1 for all i. Hence 

EE T = 9{x ) = I, 

in other words, S is orthogonal. Thus we can write 

By symmetry, ^(xm) can be diagonalizcd by an orthogonal matrix U{xm) for 
all xm- We show that U(xm) is constant. The above display yields 



x h-> 



Sp\m diag(u(x)) = *(i m )Sp\m> 
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so for all xm the eigenvectors of ^j(xm) are in the span of the columns of Sp\M j 
since rank Em ^ = p — m> rank^rcjvf). Hence the eigenvectors do not depend 
on Xm, which implies that U{xm) is constant, say equal to an orthogonal matrix 
O. Thus O t ^{xm)0 = diag(w(xA/)), for some affine vector-valued function w. 
Applying the orthogonal transformation 

^ ( T mm \ 

X ^ \ o T ) Xi 

we may assume ^(xm) is diagonal. Hence 

with E orthogonal. It is easy to show that this yields 

Ediag(^M)=( di ^^) , °™)e. 

Since an orthogonal transformation of a Brownian motion is again a Brown- 
ian motion, the E on the right-hand side can be absorbed in the underlying 
Brownian motion. The result follows. □ 

Using the above we are able to give an alternative proof of the existence and 
uniqueness results from [9], slightly strengthening the statements made there, see 
Remark 4.13 below. We first show how the well-known condition for invariance 
of the open state space M> follows from Proposition 3.7. 

Proposition 4.11. Let6(x) = diag(x) and X = R> . Then X° is stochastically 
invariant if 

Vi,Vj 7^ i : a,ij > and b{ > |. 

Proof. We apply Proposition 3.7 with $(x) = det9(x) = x\Xi . . .x p . It holds 
that 

V<f>(x)6(x) =$(x) (1 ... 1), 
and J2Li( A% T) = (! • • ■ 1 ) T > whence 

V*{x)Mx) - | =^2l[x j (a i x + b i i). 

i— 1 i j^i 

This is non- negative for all x E X = R> if > and bi > \ for all i, j. By 
applying a measure transformation with density £(\ T <r(X) -W) for some Ael p 
(which yields a probability measure by [19, Corollary A. 9]), we see that the sign 
of the diagonal elements an is irrelevant for stochastic invariance of X° . □ 
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Corollary 4.12. Suppose 9 is of the form (4-28) and assume X° ^ 0, X C 
flLii^ > 0} and {9 > 0} ^ 0. Then (2.1) is an affine SDE for a = S^H if 
and only if, up to a reparametrization of £ and v, X = n^il^i > 0} = {6 > 
0} =: T> for some q < p and 

Vi < <7,Vj < p : ftS 3 ' =0 or v t = v 3 , (4.31) 
Vi < g,Vx G <9X> 4 : /3 4 (aa; + 6) > 0. (4.32) 

Moreover, if we strengthen (4-32) to 

Vi <g,Vie &Di : ft (as + b)> i/3iX)£ T /?7, (4.33) 

i/ien I? = ni=i{ u i > 0} is stochastically invariant. 

Proof. Note that {# > 0} 7^ implies £ is non-singular. Hence 6> > if and only 
if Tr 1 e{Yr 1 ) T = diag(u) > 0, so 

v 

V=f]{v t >0}. (4.34) 

i=l 

Therefore, if X is stochastically invariant, then X = D, since 0(x) is singular 
for x G <9A? by condition (4.19), required for invariance. For the first part of the 
proposition it remains to show that conditions (4.19) and (4.20) are equivalent 
to (4.31) and (4.32). Permuting and rescaling the elements in v and column 
vectors in £ we may assume Vi = for i < q, i.e. 

X=f]{v t >0}. 

i=i 

Now condition (4.19) reads 

Vi G Q 1 j G P : (3iE J = or vj = cv t for some c > 0, 

by Proposition A. 4. The constant c can be taken equal to 1 by rescaling v and 
£. This gives the equivalence between (4.19) and (4.20), while the equivalence 
between (4.31) and (4.32) is immediate. 

For the second part of the proposition we follow the proof of Proposition 4.1 
up to Step 4. If we take 

Ui = (/3 i SS T /?7)" 1 ^, for i < q, 

then (4.10) gives 

d = B a J = u 1 (a- )- 1 (/3 l EE T /3 l T )- 1 /3 4 I]diag( V ( a;o ))I] T /37 = 1, 

by choosing xq £ T>° such that Vi(xo) = 1 for i < q. Hence the affine transfor- 
mation Lx + £ from Proposition 4.1 satisfies 

(Lx + l)i = Ui(x) = {f3 l ^ T f3j)- 1 v l (x) 7 for i G Q. 
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Let X be an affine diffusion with SDE (2.1). We have that Y := (LX + t) Q 
satisfies an affine SDE of the form 

dY t = Ji(y t )dt + dmg(^/\Yt\)dW t , 

with state space R> , where we write 

p(x) = Lp(L^ 1 (x - £)) = ax + b, 

for some matrix a with non-negative off-diagonal elements and vector b with 
non-negative components. Note that 2?° is invariant for X if and only if R> 
is invariant for Y. By Proposition 4.11, the latter holds if 6, > \ for all i, or 
equivalently, if Jii(y) > \ for all y G R> with yi = 0. Substituting y = Lx + 1 
gives 

]2i(Lx + 1) > \ for all x G dVi. 

Since pix) = Ln{L~ 1 {x - £)) and Li = (/3 i SS T /3 l T )- 1 ft, the result follows. □ 

Remark 4.13. In [9] it is assumed that q = p. Moreover, the strong existence 
and uniqueness is only proved under (4.31) and (4.33), whereas the case that 
the process X might hit the boundary dX is not treated. Note also that in [9] 
the inequality in (4.33) is strict. 

5. Quadratic state space 

In this section we consider affine diffusions where the boundary of the state 
space X is quadratic instead of linear. Let us be given a quadratic function 

$(x) = x 1 Ax + b T x + c, (5.1) 

for some symmetric non-zero A G W xp , b G W } , c G JR. We take X = {9 > 0} 
where 6 is given by (4.2) (thus X is convex) and we assume X° is a non-empty 
connected component (maximal connected subset) of {$ > 0} or {$ < 0}. 
Note that then automatically the boundary of the state space is quadratic, i.e. 
dX C {$ = 0}. By the following proposition, there are only three types possible 
for $. 

Proposition 5.1. Let X C MP be convex and assume X° is a non-empty 
connected component of > 0} or {$ < 0}, with $ given by (5.1). Then 
there exists an affine transformation such that either = x± — Y^,i=2 x i> 

= J2i=i xf + d or <&{x) = x\ — J2i=2 x i + d, for some 1 < q < p, d G R. 

Proof. Since A is symmetric, it is diagonalizable by an orthogonal matrix. By 
further scaling one can take the diagonal elements equal to —1, or 1. Using 
the equality 

x T x + b T x = {x T + \b T ){x + \b) - \b T b, 
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we can apply an affinc transformation such that for some disjoint Q, Q' C P\{1} 
the quadratic function $ is of the form 

ieQ ieQ' 



or 

ieQ ieQ' 

for some d G R. If <& is of the form (5.2), then X is of the form 

* = {a; G : a* > 5>? - £ z?}, 
ieQ ieQ' 

possibly after replacing x\ by — x\ and interchanging Q and Q' . Convexity 
of X yields that the Hessian of X^ieQ x i ~ SieQ' * s positive semi-definite, 
which implies that Q' = 0. Permuting coordinates gives Q = {2, . . . , q} with 

g = #Q + l. 

Now assume $ is of the form (5.3). We have to show that either Q' = or 
#Q < 1. Define the function / by f{x QU Q>) = E*eQ x i ~ EieQ' x l ~ d - Tncrc 
are two possible forms X can assume, namely 

X = {x G K p : xi > f{x Q uQ') 1/2 : x QU Q' G A'} 

or A" = {.r G KP : < /(x QuQ 1/2 ^QuQ' € A}, 

where if is convex with K° a non-empty connected component of {/ > 0}. In 
the first case we have that the Hessian of /(^quQ') 1 ^ 2 i s positive semi-definite, 
while in the second case the Hessian is negative semi-definite. Now suppose 
Q' 7^ 0. We show that in that case #Q < 1- Let £quQ' G K° . For i G Q' we 
have 

^! 

hence the Hessian is negative semi-definite. For i G Q we have 



#t/(*QuQ') 1/2 = -/(^QUQ')- 1/2 " z 2 /(^QUQ')- 3/2 < 0, 



^/(xquqO 172 = /(^QUQ')- 1/2 - ^/(^QuqO" 3/2 , 
which therefore also has to be negative. Now if i,j G Q and i =/= j, then 

dx,dx ] ^ XQvJQ '^'' 1 = "^^QuQ') -372 ^^' 
from which we deduce that 

d 2 d 2 



j j. I 9x 2 dxidxj 1 el \ — 1/2 

det 3 2 l 9 2 f{XQuQ>) 1 



dx,dxi dx'j 



/(^Quq0' 1/2 (/(^QUQ')- 1/2 - ^\ +^ 2 )/(^QuQ')- 3/2 ) < 0. 
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This contradicts the negative semi-definitcness of the Hessian. Thus it holds 
that #Q < 1. □ 

If (2.1) is an affine SDE with drift /x, diffusion matrix 9 and state space X = 
{9 > 0} with non-empty quadratic boundary dX C {$ = 0}, then stochastic 
invariance of X yields 

V$(x)0(x) = 0, for all x G dX. (5.4) 

by Proposition 3.1 and the remark preceding Proposition 3.4. This excludes that 
$ is of the form $(x) = Y%=i xf + d or <6(x) =x\ - J2i=2 x i+ d witn d 7^ 0. 
Indeed, suppose $(x) = + ^ with e? < (for simplicity we take q = p), 

then (5.4) reads 

x T #(x) = 0, for all x such that X^=i = 

By Lemma 5.2 below and the observation that the degree of x T £?(x) does not 
exceed the degree of $(x), it follows that 

v 

x T 9(x) = (J2 x i+ d ) cT ' for a11 x e W > 

i=l 

for some constant vector c G M p . However, since x T 9(x) has no constant terms, 
this yields x T 9(x) = for all x G MP. Lemma 5.3 below gives that 9(x) = for 
all x, which contradicts the assumption that X = {9 > 0} and dX ^ 0. Likewise 
one can show that $(x) = xf — Yli=2 x i+d with d ^ is impossible. 

Lemma 5.2. Lei $, ^> : W — >• K be multivariate polynomials and O C MP a 
non-empty open set such that for x G O we have ^(x) = whenever $(x) = 0. 
Suppose in addition that for x G O, $(x) is of the form 

d 

$(x) = JJ^i -rj(x P \ { i } )), 

3=1 

/or some e? G N and some Tj{xp\si\) G M which are mutually different for all 
x G O. TTiera ^(x) = C(x)$(x) /or some multivariate polynomial C . 

Proof. Viewing x\ i— > $(x) as a univariate polynomial in xi with distinct roots 
r i( :£ P\{i}) f° r nxe d ip\{i}, it follows by the factor theorem for polynomials that 

\&(x) = C(x)<&(x), for x G O, 

for some function C that is a polynomial in x\. Since \1/ is a multivariate poly- 
nomial that only depends on x\ by the term C(x)xf, it follows that C is a 
multivariate polynomial. The equality in the above display can be extended to 
W by applying a uniqueness theorem for holomorphic functions, see [18, p. 226]. 
This yields the result. □ 
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Lemma 5.3. Let 8 : W -> W xp : x i-> A" + J2l =1 A k x k with A k e «p x p 
symmetric and assume x T 8(x) = for all x. Then 8(x) = for all x. 

Proof. It is clear that A = 0. We show that A k = for k > 0. It holds that 

p p p 

o = (x T ^ k x k ). = J2J2 x ^a% = J2 x i x M% + A kj) 

k=l k=l i=l l<i<k<p 

i=l 

for all j. Hence for all i, j we have A], = and A k , = — A k - for k ^ i. Since A fc 



is symmetric, the latter gives A k , = —A\, = — A*- fc . So if we permute the indices 



ij ij kj 

— A i — — A i 
Hj ~ A kj — A jk- 

k by the cycle (i i— > j, j n- fc, fc >— > i), then A^- gets a minus sign. Permuting 
the indices repeatedly we obtain 



A k — — A 1 — A 3 — — A k 
A ij ~ A jk — A ki — A iv 

which implies A k j = for all i,j and k ^ i. Hence A k = for all k, as we have 
already shown that A\j = for all i, j. □ 

Thus in order to characterize all affine diffusions with quadratic state space, 
there are two cases to consider, namely &(x) = x\ — J2i=2 x i anc ^ = 
x\ — ^21=2 x f ■ ^ n ^ ne nex t subsections we characterize for these two forms of <E> 
all possible 8 which can act as a diffusion matrix of an affine SDE, i.e. which 
8 satisfy (5.4). Moreover, we are able to construct a square root a of 8 such 
that (2.1) is an affine SDE with quadratic state space X, generalizing the 2- 
dimensional setting as treated in [8, Section 12] and [12]. In particular we show 
existence and uniqueness of a strong solution. 

5.1. Parabolic state space 

Assume $ is of the form $(cc) = x\ — J2i=2 x h with 1 < q < p. The state space 
X then necessarily equals X = > 0}. For x € W we write x = (xi,y, z) £ 
R 1 x R 9_1 x M. p ~ q and we define affine matrix- valued functions £ and r\ by 

t(x)=( 4Xl 2yT ) n(x)=( ° ° ° 

UJ \2y I )' 1[ ) \T 12 (y) T 13 (y) ... r,_ a , g _i(y) 

with Tij : R"- 1 -> for 1 < i < j < q given by T <3 -(y)< - (</),■ = 

Tij{y)k = for fc ^ i, j. Moreover, we write Q = {2, . . . , g}. Now, condition (5.4) 
reads 

(l -2y T 0) 8(x) = for all x e R*> such that an = y T y. (5.5) 

We use the following lemmas to characterize those 8 that satisfy this condition 
in Proposition 5.6. 
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Lemma 5.4. Consider the linear space 

C = {a : R p -> R q affine | (l -2y T ) a(x) = for all x with x-y = y T y} . 

(5.6) 

Then a basis for C is formed by the columns of C and r\. 

Proof. Clearly these columns are linearly independent elements of C. To prove 
that they span C we use a dimension argument. Let Aff(R p ,R l? ) denote the 
space of affine functions from W to R q and let Quadr(R p , R) / {x\ — y T y) be the 
space of quadratic functions from R p to R, modulo x\ — y T y (that is, p and q 
arc equivalent if p(x) — q(x) = c{x\ — y T y) for some constant c). Consider the 
linear operator 

L : Aff(R p , R 9 ) ^ Quadr(R p ,R)/(x 1 - y T y) : a(x) i— > (l -2y T ) a(x), 

and note that C = kerL, in view of Lemma 5.2. By the dimension theorem for 
linear operators, we have 

dim Aff(R p , R 9 ) = dim ker L + dim imL. 

It holds that dim Aff(R p , R f/ ) = pq + q. Since x\ = y T y, a basis for im L is given 
by 

{l,x 2 , ■ ■ ■ ,x p } U {xiX 3 : 2 < i < q, 1 < j < p}, 

whence 

'q-A _ _ (q-l 



dim imL = p+ (q - l)p - I ^ j = pq . 2 

It follows that dimkerL = q+ |(g — l)(g — 2), which is the number of columns 
in C and rj. Thus the columns span the kernel of L. □ 

Lemma 5.5. Let C be defined by (5.6) and suppose M : W — s- R qxq is affine 
and M(x) is symmetric for all x € R p . // the columns of M are in C, then 
AI = c£ for some cel. 

Proof. By Lemma 5.4 there exist matrices A and B such that 

M(x) = £(x)A + r)(x)B. 

Write T(y) = (L ij (y)) 1 < i<J<g and B = (p, 1 Bj. Then the above display reads 



M{x) 



( Ax x A lx + 2y T A Q1 A Xl A 1Q + 2y T A QQ _ 

\2yA 11 + Aqi + T(y)B 1 2yA 1Q + A QQ + T(y)B 



Since M(x) is symmetric it immediately follows that Aiq = and Aqi — 0. 
Define N(x) = M(x) — An<^(x). Then N is also symmetric. We have 



N(x) = 



2y T C 
T{y)B 1 C + T(y)Bj' 
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with C = A QQ - AnI. This yields C = C T and T(y)B 1 = 2Cy. Since y T T(y) = 
0, the latter implies y T Cy = 0, whence C = 0, as C is diagonalizable by an 
orthogonal matrix. Thus Aqq = A\il and it remains to show that B = 0. 

It holds that T(y)B is symmetric and y T T(y)B = 0. Lemma 5.3 yields 
T(y)B = 0, whence B = by linear independence of the columns of T(y), 
as we needed to prove. □ 

Proposition 5.6. If (5.5) holds, then necessarily 9 is of the form 

•M-C$# &))■ <"> 

/or some c > 0, tw't/i -A(x) = Q{x)A\ + r]{x)A2 for some matrices A\^A<2 and 
B :W -> m(p-?) x (p-9) a^ne and symmetric. Moreover, if c = 1 ant! A x = 0, 
i/ien ii ZioWs </ia< 

B(s) - Ajr ? (a;) T 77(a;)A 2 > 0, (5.8) 

/or aZZa; S X = {9 > 0}. 

Proof. The first part follows from Lemma 5.4 and Lemma 5.5. It remains to 
show (5.8). 

Suppose c = 1 and A\ = 0. By positive semi-definitcness of 9, we have 

< (v T w T ) 9{x) 

= v T C,(x)v + 2v T n(x)A2W + w 1 B(x)w, 

for all v G R 9 , w G R p ~ 9 , sc G <Y. Fix to G R^ 9 , a: G X arbitrarily and take 
v = —rj(x)A2W. Noting that £(21)77(2;) = n(x) for all x G R p , the above display 
then reads 

w T B(x)w — w T Aj n(x) T r)(x)A2W > 0, 
which proves (5.8). □ 




To show the existence of an affinc diffusion with parabolic state space in 
Theorem 5.8, we need the following result. Its proof is based on a modification 
of a result by Yamada and Watanabe [21, Theorem 1]. 

Proposition 5.7. There exists a unique strong solution to the SDE 



Xi,t\ _ ( a\iX\,t + a 1Q Yt + h 




1 y,) = { Z Q y,7l' Q 

Proof. By continuity of the coefficients and satisfaction of the linear growth 
condition (2.2), there exists a weak solution {{X\, Y), W) on some filtered prob- 
ability space (O, J 7 , (Jt),P) carrying a Brownian motion W. For strong exis- 
tence and uniqueness it suffices to prove pathwise uniqueness. Therefore, assume 
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((X\,Y), W) is another solution on the same probability space. We see that the 
equation for Y does not contain X\ and in fact it is an SDE which has a unique 
strong solution, whence Y t — Y t a.s. for all t > 0. Write Z = X\ — Y T Y and 
Z = X 1 - Y T Y. Then it follows that 

d(Z t - Z t ) = a n (Z h t - Z ht )dt + 2 (y\z7\ - \f\Z^J dW ht . 

Arguing as in the proof of [15, Proposition 5.2.13] (for instance), we deduce that 
Z t = Z t a.s. for all t > 0. Thus X x ,t = Z t + Y^Y t = Z t + Y t T Y t = X ht a.s. for 
all t > 0. □ 

Theorem 5.8. Let $(x) = x\ — y T y and suppose {6 > 0} = {$ > 0}. Then 
there exists an affine SDE with drift \i, diffusion matrix 6 and state space X = 
{9 > 0} if and only if 

V$(x)6(x) = 0, (5.9) 
V$(a)/i(x) >tr(e QQ (x)) (5.10) 

for all x G dX = {$ = 0}. 

Proof. Conditions (3.6) and (3.7) reduce to (5.9) and (5.10). The "only if" -part 
follows, as these boundary conditions are necessary for stochastic invariance of 
X by Proposition 3.4. 

Suppose the boundary conditions (5.9) and (5.10) hold. Then by Proposi- 
tion 5.6 it holds that 6 is of the form (5.7) for some c > 0. If c = 0, then 
necessarily X = {B > 0} and we take 



, v _ (0 
a[X) ~ \0 \B(x)\^ 2 



It is easy to see that (2.1) admits a unique strong solution and that X is invari- 
ant. If c > 0, then we apply the linear transformation x\ i— > c~ 1 x±, y i— > c _1//2 y, 
z i y z c~ l Ai(x\, y), so that we may assume c = 1 and A\ = 0. Then (5.8) 
holds and as a square root of 9 on X we take 



a{x) 

with £ and p defined by 



m o 

Aj v (x) T p{x) 



P=\B-AJ V T V A 2 \^ 2 . 

To see that a(x)a(x) T = 9(x) for x E X, note that £(x)r)(x) = rj(x) for all 
x G A". Proposition 3.4 gives that A is invariant for (2.1). It remains to prove 
existence and uniqueness of a strong solution. 
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The boundary condition (5.10) reads 

fL\{x) — 2y T fig(y) > q — 1, for all x such that x\ = y T y, 
so necessarily fi{i}uQ admits the form 




(5.11) 



Proposition 5.7 gives the result. 



□ 



For stochastic invariance of a parabolic state space X, we now give sufficient 
conditions on the diffusion matrix and sufficient and necessary conditions on 
the drift, analogous to the admissibility conditions for polyhedral state spaces 
in canonical form. 

Proposition 5.9. Condition (5.8) holds if B admits the form 



The right-hand side equals r](x) J rj(x) , while the left-hand side is equal to (q — 
2) (y2 i yf) I. which is smaller than (q — 2)x±l for x € X. Hence 



(q - 2) Xl A 2 f A 2 - Al r](x) 1 r](x)A2 = A± ((q - 2) Xl l - rj(x) 1 rj(x))A 2 > 0, 



Proposition 5.10. Suppose X is an affine diffusion with state space X = {9 > 
0} = {x £ M. p : x\ > y T y} and 8qq ^ 0. Then X is a linear transformation of 
an affine diffusion with the same state space X , diffusion matrix of the form 
(5.7) with c — 1 and with drift fi(x) = ax + b satisfying (5.11) as well as 



B{x) = (q-2) Xl AjA 2 . 



Proof. One can show using Cauchy-Schwarz that 



7](x)i7](x)i I > rj{x) i T)(x)i, for all x and i. 



Hence 




for x £ X, which yields the result. 



□ 



anl - 2a QQ = diag(d Q ) 



(5.12) 
(5.13) 
(5.14) 
(5.15) 

(5.16) 



d Ql > 
d Q2 = 
a\Q 2 = 2bQ 2 




for some vector d and some disjoint Q\ and Q 2 with Q = Q\ U Q 2 
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Proof. As in the proof of Theorem 5.8 we can assume 9 is of the form (5.7) with 
c = 1. Then tr (0qq(x)) = q — 1 for all x, so the boundary condition (5.10) for 
the drift reads 

y T (anI - 2a QQ )y + (a 1Q - 2b T Q )y + h - q + 1 > 0, for all y e W 1 . 

For this it is necessary that M := anl— 2agg is positive semi-definite. Moreover, 
if y is in the kernel of M, then y should also be in the kernel of a%Q — 2bg. We 
can diagonalize M by an orthogonal matrix O, so D = OAIO T is diagonal with 
positive diagonal elements di for i € Qi and di = for i € Q 2 = Q\Q\, for some 
Qi C Q. Applying the orthogonal transformation y t— > Oy, the above condition 
becomes 

dixj + ^2 (on - 2hi)xi + bi - q + 1 > 0, for all x Ql . 

ieQi ieQi 

We can write the left-hand side as 

53 Mm + \dT l { ail ~2bi)) 2 - JJ] t^^Oii - 2bif + h - q + 1, 
ieQi ieQi 

which is non-negative for all xq ± if and only if 

— \ J2 C X («H - 2bi) 2 + h-q + l>0. 

This yields the result. □ 

Remark 5.11. Note that (5.12) implies that oqq is diagonal. Hence the coordi- 
nates of Xq are mutually independent. 

Proposition 5.12. Consider the situation of Proposition 5.10. If we strengthen 
condition (5.16) to 

bi >q + l+ ]T \dr 1 (a ll -2b l ) 2 , (5.16') 
ieQi 

then X° becomes invariant. 

Proof. It suffices to verify condition (3.9) of Proposition 3.7. Recall that $(x) = 
x\ — y T y and 6(x) = A + Y^h=i A x Xi is of the form (5.7) with c = 1. It follows 
that V^(x)(n(x) - ^Lii^Y) c q uals 

a n xi + a 1Q y - 2y T a QQ y - 2y T b Q + 6i - q - 1. 

If in addition to conditions (5.12) - (5.16') also an > is imposed, then for 
x G X the above display is bounded from below by 

y T (a n I - 2a QQ )y + (a 1Q - 2bg)y + b 1 -q-l. 
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This is non-negative for all y G R 9 under the imposed assumptions, similar as 
in the proof of Proposition 5.10, which yields (3.9). The non-negativity of an 
can be dispensed with, as shown as follows. 

By applying a measure transformation with density £(X T a(X) ■ W) for some 
A G M p with A; = for i ^ 1 (which yields a probability measure by [19, 
Corollary A.9]), we see that X is also invariant for the SDE with drift p(x) = 
ax + b + 9(x)X = ax + b, where on = an + 4A, clqq = oqq + 2AI and the 
remaining coordinates unaltered. Note that an I — 2a qq = an I — 2oqq, so 
conditions (5.12) - (5.16') are not affected by such a measure transformation. 
This gives the result. □ 



5.2. Conical state space 

Let p > q > 1. For x G MP wc write x = (xi,y,z) G K 1 x W'- 1 x W~ q . We 
consider the quadratic form 

where Q = {2, . . . , q} and define affine matrix-valued functions £ and r\ by 

C(a ° = (y xii)' v{x) = (r 12 °(y) r 13 °(y) r,_ a °_ 1 (i,)) ' 

with T y : R*- 1 -> K*- 1 for 1 < * < j < q given by T l3 {y) % = Vj , T ij {y) j = -y u 
Tij(y)k = for k ^ By applying a reflection, we may assume the state 
space X is of the form X = {3> > 0} n {x\ > 0}. Analogously to Lemma 5.4 and 
Lemma 5.5 we have the following. 

Lemma 5.13. Consider the linear space 

C = {a : W ->• R 9 affine \ (x t -y T ) a(x) = for all x with x\ = y T y} . 

(5.17) 

Then a basis for C is formed by the columns of £ and r\. 

Proof. Similar to the proof of Lemma 5.4. □ 
Lemma 5.14. Consider the linear space 

M = {M : M p R qxq affine | M (x) symmetric and AP G C for all x, i] 
with C defined by (5.17). Then a basis for Ai is given by 

B = {(,p(l),...,p(q-1)}, 
with p(i) an affine symmetric-matrix valued function defined by 

p(i) i+ i y T ) , 

p(i)n :x^yi, 

p(i)jj : x h-> -j/,, for j ^ l,i + 1, 
p(i)jk ■ x 0, if j,k^i and j ^ k. 



P. Spreij and E. Veerman/Affine diffusions with non-canonical state space 



32 



Proof. Clearly the elements of B are linearly independent elements of A4. It 
remains to show that they span A4. 

Let MeMbe arbitrary. By Lemma 5.13 there exist matrices A and B such 
that 

M(x) = £(x)A + T)(x)B. 



M(x) 



Write Q = {2,...,q}, T{y) = (T (j -(y))i< i<J<g and B = {^B 1 B). Then the 
above display reads 

x 1 A 11 +y T A Q1 x 1 A 1Q +y T A QQ \ 

K yA xl + Xl A Ql + T(y)B 1 yA 1Q + Xl A QQ + T(y)Bj 

Symmetry of M(x) yields 

Aiq = A Q1 , 

A QQ = A QQ' 

yA\i + T(y)B 1 = Ag Q y : 

yA 1Q +T(y)B = (yA 1Q +T(y)B) T . 

Since y T T{y) = 0, the second equation together with the third gives 

= y T T(y)B 1 = y T (A QQ - A lx l)y, 

which implies Aqq — Axil = 0, as Aqq — Ai±I is symmetric and thus diagonal- 
izable by an orthogonal matrix. Define 

N = M-A 11 (-J2 a uP(})- 

Then N <G M and N is of the form 

N(x) fQ ° 



EkeQ&VkJ' 

for some symmetric ((q — 1) x (q — l))-matrices C k . By Lemma 5.3 it follows 
that TV = 0. □ 



Unlike the parabolic case, for a general closed conical state space we are not 
able to find a square root such that strong existence and uniqueness for the 
resulting SDE can be proved. An exception is the two-dimensional cone, as this 
is just a polyhedron which has already been covered in Section 4. The following 
example shows that problems appear for closed cones in higher dimensions. 

Example 5.15. For p = q = 3, a basis for M. is given by 

fxi yi y 2 \ 

C(X) =\ yi X, ) , p(l)(x) : 

P (2)(X) : 



Xl 


in 


.'72 


yi 


■>'i 





!)2 








tji 


Xl 
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Xl 
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Xl 


!)i 
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Note that not only ( but also £ + p(l) and ( + p(2) arc positive semi-definite on 
X = {$ > 0} n {.ti > 0} = {x € M 3 : xi > 0, x\ > y T y}. The structure of these 
matrices appears to be too complex to compute a manageable square root. 

However, Proposition 3.7 enables us to derive sufficient conditions for stochas- 
tic invariance of the open conical state space {£ > 0}. This can be used to show 
existence of a unique strong solution for the afhnc SDE (2.1) with square root 
a = ICI 1 / 2 , see the next proposition. Note that this approach is not applicable 
for ( + p(l) and £ + p(2) in Example 5.15, as these matrices are singular on 
the whole of R 3 . We leave the question of existence of an afhne diffusion with a 
closed conical state space open for further research. 

Theorem 5.16. There exists an affine SDE with drift fi(x) = ax + b, diffusion 
matrix 



existence and uniqueness for (2.1) follows (see [15, Theorem 5.2.5]) as soon as we 
have shown stochastic invariance of X. It holds that X is a connected component 
of {<£> > 0}, with $(x) = x\ — y T y. Therefore, in view of Proposition 3.7, 
it suffices to prove (3.9) and (3.10). The first condition is immediate. For the 
second one, a calculation shows that VQ(x)(fi(x) — \ X^i^*) 1 ) equals 

2{a xl x\ + x 1 (a 1Q - Og^y + (bi - \p)x\ - y T a QQ y - b^y). 

This is non-negative for all x G X if (5.18) - (5.20) hold and an > 0. Indeed, in 
that case we have for x e X = {x\ > y T y} n {x\ > 0} that 



since — (6g,y) > — ||6q||||j/|| > — ||6q||xi by Cauchy-Schwarz. The non-negativity 
of an can be dispensed with, by the same arguments as in the proof of Propo- 




and conical state space X = {9 > 0} = {x\ > (y T y) 1 ^ 2 } if 




anx\ + X!(a 1Q - Og-Jy + (h - \p)x t - y T a QQ y - b^y 

> y T (anI - a QQ )y + (6 X - \p)x^ - (b Q , y) 

> y T (anI - a QQ )y + (bi - \<p - \\b Q \\)x 1 > 0, 



sition 5.12. 



□ 



Appendix A: Convex analysis 



In this section we state and prove the results on convex analysis applied in 
Section 4. Let X be given by (4.1) and in addition to an affine function u we 
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are given an affine function d by 

d : R p -> R : x i-> ax + b, 

for some a € R lxp , b G R. Proposition A.l below is the main result, which yields 
Proposition A. 2 to tackle the drift and Proposition A. 4 to tackle the diffusion 
matrix of afhnc diffusions with non-canonical polyhedral state space. 

Proposition A.l. Suppose X C {d > 0}. Then there exist c > and A G R>o 9 ; 
such that 

d = Xu + c. 

Proof. We give a proof by contradiction. Let 

K = { (A 7 , X8 + c) : X G M> x 9 , c > 0}. 

Suppose (a, 6) G" /C. Since IC is a closed convex set, (a, 6) is strictly separated 
from K, by the Separating Hypcrplane Theorem. Therefore, there exist j/el p 
and i/o€l such that ((y, y ), (k, ho)) > ((y, yo), (a, b)} for all (fc, fco) G K,, i.e. 

fcy + koyo > ay + byo for all (fc, fco) G /C. 

In other words, for all A; > and c > we have 

^2 K{liV + hyo) + cy > ay + by . 

i 

It easily follows that 

ay + by o <0 (A.l) 
HV + Siyo > (A.2) 
yo > 0. (A.3) 

Using this we construct x G X for which d(x) < 0. Suppose yo > 0. Then we 
take x = y/yo- Indeed, Ui(x) = (jiy + Siyo) / yo > 0, so x £ X. But d{x) = (ay + 
byo)/yo < 0, which is a contradiction. Suppose yo = 0. Then we take an arbitrary 
xo G X and let xn — xo + Ny, with N G N. Then Ui(xn) = Ui(xo) + N"/iy > 
for all N, so x^ G X, but d(x^) = d(xo) + Nay < for N big enough. □ 

Proposition A.2. Suppose dXi C {<i > 0} for some i E Q. Then there exist 
c > and X G R lx? iwi/i A.,- > /or j G suc/l 

d = Au + c. 

Proof. Let it := — Uj. Then dA^ = Q* =0 {uj > 0} and d(x) > for x G 
Hence we can apply Proposition A.l, which gives the existence of A^ > with 
j =0,...,q and c > such that 

with Aj = Xj > for j ^ i and A; = Xi — Xo- □ 
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Lemma A. 3. Assume Q is minimal. It holds that dXi ^ for all i £ Q. 

Proof. Fix i < q. By minimality of Q we can choose x £ W such that Ui(x) < 
and > for all j ^ i. Since X ^ 0, we can choose y <G X. Then > 

for all j. For t £ [0, 1] it holds that 

+ (1 - = *Wj(a;) + (1 - t)uj(y), 

which is non-negative for j =/= i. For t = Ui(y) / (ui(y) — Ui(x)) we have Ui(tx + 
(1 - t)y) = 0, so tx+(l- t)y £ 9AT f . □ 

Proposition A. 4. Assume Q is minimal. Suppose dXi d {d = 0} for some 
i £ Q. Then there exists Xi £ M. such that v{x) = XiUi(x) for x £ X. If X° =/= 0, 
then v(x) = XiUi(x) for all x £ M. p . 

Proof. We have dXi C {v > 0} and dXi C { — v > 0}. Applying Proposition A. 2 
with d = v respectively d = — v, we derive that 

9 

v(x) = XjUj(x) + c\ 

3 = 1 

i 

—v(x) = }]fj,jUj(x) +c 2 , 
i=i 

for some A,/i € R lxp with Xj,fij > for j =^ i and Ci,C2 > 0. Adding the 
equations in the above display gives 

9 

= y^(Aj + fXj)uj(x) + c\ + c 2 . 
j'=i 

By Lemma A. 3 we can choose x £ dXi and deduce that c\ = c 2 = 0. So 
-(A, + Hi)v,i(x) = y^(Aj + fj,j)uj(x). 

By minimality of Q we can choose x € K p such that «i(x) < and Uj(x) > 
for all j i. This gives that c := A^ + \Xi > 0. If c > 0, then for a; £ we have 

< = -c- 1 £)(Aj + Mi)%0) < 0, 

whence Ui(x) = for x £ So A" = dXi C = 0} and v{x) — Ui(x) = for 
x £ X. If c = 0, then ^-^(Aj + fj,j)uj(x) = for all x. This holds in particular 
for x £ i.e. for a; such that Uj(x) > for all j. Hence for x £ X we have 
XjUj(x) = jjtjUj{x) = for all j 7^ i, so 

9 

v(x) = ^^Aj-u^x) + cj = Ajitj(x), (A. 4) 

for a; S Af. If Af° 7^ 0, then choosing x <G Af° gives %(x) > for all j, which 
implies X 3 = for all j i. Then (A.4) holds for all x £ W. □ 
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